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ABSTRACT 

We present interferometric observations resolving the CO emission in the four gas-rich lenticular 
galaxies NGC 3032, NGC 4150, NGC 4459, and NGC 4526, and we compare the CO distribution and 
kinematics to those of the stars and ionized gas. Counterrotation documents an external origin for the 
gas in at least one case (NGC 3032), and the comparisons to stellar and ionized gas substructures in 
all four galaxies offer insights into their formation histories. The molecular gas is found in kpc-scale 
disks with mostly regular kinematics and average surface densities of 100 to 200 Mq pc~^. The disks 
are well aligned with the stellar photometric and kinematic axes. In the two more luminous Virgo 
Cluster members NGC 4459 and NGC 4526 the molecular gas shows excellent agreement with circular 
velocities derived independently from detailed modeling of stellar kinematic data. There are also two 
puzzling instances of disagreements between stellar kinematics and gas kinematics on sub-kpc scales. 
In the inner arcseconds of NGC 3032 the CO velocities are significantly lower than the inferred circular 
velocities, and the reasons may possibly be related to the external origin of the gas but are not well 
understood. In addition, the very young population of stars in the core of NGC 4150 appears to have 
the opposite sense of rotation from the molecular gas. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: ISM — galaxies: kinematics 
and dynamics — galaxies: structure — galaxies: evolution — galaxies: individual 
(NGC 4459, NGC 4150, NGC 3032, NGC 4526) 
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1. INTRODUCTION 

Early-type galaxies, the ellipticals and lent iculars, are 
gener ally poorer in cold gas than spirals (e.g. iLees et all 
Il99l . This fact is, of course, ultimately responsible for 
the color difference between early- and late-type galax- 
ies and their locat ions in the red sequence or the blue 
cloud (e.g. Bald r v et al.ll2004f ). However, cold atomic 
and molecular gas are not entirely absent from all early- 
type galaxies, and therein lie some important clues to 
both the past and the future of early-type galaxies. The 
origin of the cold gas in early-type galaxies can serve as 
a tracer of their assembly histories, and the properties 
of the cold gas (the raw material for star formation) of- 
fer insights into possible morphological and dynamical 
e volution through sta r form ation. 

iFaber fc Gallagheil (|1976( ) predicted that mass loss 
from evolved stars in early-type galaxies should pro- 
duce detectable amou nts of gas over a Hubble time. 
More recent work (e.g. ICiotti et al.lll991h makes predic- 
tions which are different quantitatively but not qual- 
itatively. Of course, it is likely that the gas recy- 
cled into the interstellar medium will have a complex 
thermal history, including heating to X-ray tempera- 
tures and possible cooling all the way to the forma- 
tion of molecules (jBrighenti k Mathewd I1996L Il997t ). 
But since the stellar mass l oss i s unavoidable and 
iTemi. Brightenti. fc Mathews! ()2007f ) have shown that 
the dusty ejecta can cool to molecular temperatures in 
galaxy centers, it is very natural to expect that some 
early-type galaxies could contain cold gas. In fact, as 
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iTemi. Brightenti. fc Mathews! ()2007f ) have alluded, the 
question of why there are early-type galaxies without cold 
gas may be just as interesting as why there are early-type 
galaxies with cold gas. 

However, if the cold gas ultimately originated in 
the stars, one might expect a correlation between the 
luminosity of an early-type galaxy and its cold gas 
content, which has not be en seen (Wardle fc Knao: ' 



19861 iHuchtmeier Sage fc He nkel 1995; Lees ct al. 199^ 
Knapp fc RupenI 119961: ICombes'et aLH2007t) . The lack 
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of such a correlation is then often used to suggest that 
whatever molecular gas does exist in early-type galax- 
ies is completely unrelated to the internal stellar mass 
loss, having been acquired from external sources such 
as a satellite galaxy or perhaps even the intergalactic 
medium. In a slightly more sophisticated analysis of 
both the atomic and mo l ecula r gas contents of lenticular 
galaxies, ISage fc WelchI (|200(j ) hypothesized a variation 
on this theme in which the atomic gas could be of exter- 
nal origin whereas the molecular gas could be of internal 
origin. 

Due to the complex nature of the thermal evolution 
of the stellar mass loss, plus environmentally dependent 
interaction with an intracluster medium and potential 
time-dependent effects of AGN feedback on the ISM, we 
argue that the total gas masses by themselves do not 
offer compelling answers about the origin of the cold 
gas in early-type galaxies. However, a comparison of 
the gaseous and stellar angular momenta should provide 
much stronger constraints. The distribution and kine- 
matics of the gas can reveal recent gravitational inter- 
actions, and gas which is counterrotating with respect 
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to the stars almost certainly did not originate in those 
stars. Maps which resolve the molecular gas are thus 
crucial tools which help us to read the assembly histories 
of early- type galaxies. 

Molecular gas in early-type galaxies is also particularly 
interesting because such gas is the raw material for star 
formation activity. In recent years there have been sug- 
gestions that current-day star formation is takin g place in 
as many as 30% of early-t ype galaxies at z = (|Yi et al] 
l2005t iKavirai et al.ll2006[ ). Such inferences are at least 
superficia lly consistent with t h e CO detection rates 
found by iSage fc WelchI j2n0^ . iSage. Welch, fc Yound 
(|2007| ). and lCombes et al.l ([2003). But the distribution of 
the molecular gas will determine whether the star forma- 
tion takes place in a nuclear starburst or in an extended 
disk and, therefore, it will also determine the efficiency 
with which AGN feedback could disrupt the star forma- 
tion. Testing our theoretical understanding of star for- 
mation in early-type galaxies also requires estimates of 
the molecular disk sizes and gas surface densities. Thus, 
resolved maps of molecular gas hold the keys to under- 
standing the morphological evolution of early-type galax- 
ies through star formation. 

We begin to address these issues on the past and fu- 
ture of early-type galaxies through interferometric maps 
which resolve the CO emission in four nearby lenticular 
galaxies. Total molecular masses are comm only available 
for SOs, but resolved C O maps are rare fe.g. lOkuda et al.l 
l2005t iDas et aLll2005D and the four presented here repre- 
sent a significant addition to the literature on the subject. 
In all four cases the molecular gas is found in kpc-scale 
disks, and in one case the molecular gas counterrotates 
with respect to the stars. We also present comparisons 
between the CO kinematics and circular velocity curves 
derived independently from stellar kinematic data. We 
discuss the origin of the gas and its implications for mass 
modeling. A forthcoming paper will also present compar- 
isons to the stellar populations in the galaxies and their 
ionized gas properties, with special emphasis on the dis- 
tribution and kinematics of the young stellar population, 
showing clear evidence for disk growth through star for- 
mation. 

2. SELECTION AND PROPERTIES OF THE GALAXIES 

The four galaxies NGC 3032, NGC 4150, NGC 4459, 
and NGC 4526 were observed because of their 
membership in the SA URON integral-fi e ld su rvey of 
early- type galaxies (see Ide Zeeuw et al.l l2002l ). which 
means that detailed maps of the stellar kinemat- 
ics, ionized-gas kinematics, and stellar populations 
are available for the regions within an effective 
radius (e.g. Qj lmscllcm et al] |2004 fSarzi et all 1200^ : 
iKuntschner et al. 2006i'). They also ha,ye CO detections 
(IThronson et al. [119891: ISage fc WrobellHosl llees" et al.l 
ll99lHCombes' et al making them suitable for high 

spatial resolution comparisons of the stellar and molec- 
ular disks. Basic parameters of the galaxies are summa- 
rized in Table [H 

NGC 3032 is a relatively faint {Mb = -18.8) field 
lenticular with complex stellar kinematics. Having high 
angular momentum per unit mass, it belongs to the so- 
called fast rotator class (jEmsellem et al. | |2007f). and has 
a small counter-rotating core (R » 2' MMcDermid et al.l 
|20063). The ionised gas is co-spatial with dusty spiral 



arms extending to ab out an effective r adius Re and is 
also counter-rotating (|Sarzi et al.l I2005D . H/J emission 
is the strongest in the SAURON sample while the ratio 
[O Illj/H/J is weakest (except in the very center). The 
absorption linestrength distributions are unusual, sug- 
gesting young stars everywhere with a lower metallic- 
ity a nd a-enhancement in the center (jKuntschner et all 
120061: Kuntschner et al. 08, in prep.), probably indicating 
a recent starburst. 

NGC 4150 belongs to the Coma I Cloud and is 
low-luminosity {Mb — —18.5) lenticular in many re- 
spects similar to NG C 3032. It is a fast rotator 
(jEmsellem et al.l l2007t ) with moderate rotation and a 
dip i n the central velocity dispersion (jEmsellem et al.l 
'2004!) , and also harbours a central counter-rotating core 
(.McDermid et all l2006bD where messy dust is found. 
Contrary to NGC 3032, however, the extended ionised 
gas is co-rotating with the bulk of the stars, although 
slightly misaligned (jSarzi et al.ll2005l ) . The linestrengths 
again indicate young stars everywhere, particularly in the 
center, but the metallici ty and a elements are rather flat 
(jKuntschner et aLll2006t Kuntschner et al. 08, in prep.). 

NGC 4459 is an average {Mb = -20.0) Virgo lenticu- 
lar with rapid rotation and the pinch of the stellar isove- 
locity contours characteristic of decoupled central disks 
(Emsellcm et al. 2004). This disk is co-spatial with a reg- 
ular but flocculent dust disk and rapidly rotating ionised- 
gas with low [O IIIj/ H/? ratio (except in the very cen- 
ter; jS^iTitlaLjllOOH)- The Hnestrengths indicate much 
younger stars and higher metallicity in the central disk 
(Kuntschner et al. 2006; Kuntschner et al. 08, in prep.). 

NGC 4526 is a Virgo lenticular similar to NGC 4459 
but slightly brighter {Mb = — 20.7). It harbours what 
is probably the strongest decoupled central stellar disk 
in the SAURON sample (observed nearly edge-on ), em- 
bedd ed in a more slowly rotating bulge (EmscUem et al.l 
l2004f ). The bulge has a slight triaxiality which suggests 
the presence of a bar. A rather large but regular dust and 
ionised-gas disk is co-spatial with th e stellar disk, an d 
has uniformly low [O III]/H/3 ratio (jSarzi et al.ll2005D . 
The linestrengths clearly show that the central disk is 
much younger and metal rich than the bulk of the galaxy 
(jKuntschner et aLll2006 : Kuntschner et al. 08, in prep.). 

The four galaxies studied here thus fall naturally 
into two classes. NGC 3032 and NGC 4150 are rela- 
tively faint and have low velocity dispersion, but har- 
bour small counter-rotating cores and have extended 
marginally ordered gaseous disks with pervasive young 
stars. NGC 4459 and NGC 4526 have rapid rotation 
but even more rapidly co-rotating decoupled (and spa- 
tially well-constrained) central stellar disks, co-spatial 
with regular dust disk populated with young stars. The 
former thus point to a more troubled recent past, per- 
haps involving accretion or minor mergers, while the lat- 
ter suggest a more peaceful immediate past dominated 
by secular evolution. The current molecular gas obser- 
vations should thus allow to build on these data and test 
those ideas further. 

3. CO DATA 

NGC 3032, NGC 4150, NGC 4459, and NGC 4526 were 
observed in the ^^CO J = 1 — > line with the 10-element 
Berkeley-IUinois-Maryland Associati on (BIMA) rnillime - 
ter interferometer at Hat Creek, CA (jWelch et al.lll99(j) . 
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TABLE 1 

Sample Galaxies - Basic Properties 







Nrr* zLi 






RA (J2000.0) 


09 52 08.2 


12 10 33.6 


12 29 00.0 


12 34 03.0 


Dec 


+29 14 10 


+30 24 06 


13 58 43 


07 41 57 


Velocity (km s^^) 


1533 (5) 


226 (22) 


1210 (16) 


575 (24) 


Distance (Mpc) 


21.4 


13.4 


15.7 


16.4 


Type 


SABOO(r) 


SOO(r)? 


SO+(r) 


SABOO(s): 


T 


-1.7 


-2.4 


-2.0 


-1.6 


Mb 


-18.8 


-18.5 


-20.0 


-20.7 


{B-V)e 


0.63 


0.83 


0.97 


0.98 


Rc (") 


9 


18 


35 


44 


CTQ (km s ^) 


82 


148 


174 


256 


CO flux (Jykms-i) 


93 (18) 


26 (5) 


56 (11) 


180 (36) 


H2 mass (10® M0) 


5.0 (1.0) 


0.55 (0.11) 


1.6 (0.3) 


5.7 (1.1) 



Note. — Data are taken from NASA's Extragalactic Database (NED) and the Lyon Extragalactic Database 
(LEDA). The velocities are stellar m easurements with the exceptio n of NGC 3032's, which is a HI velocity. Distanc e 
measurements are from [Tbnry et all (20011) and IMei et al.l II2005I) and Mb are taken from [Emsellem et al.l II2007I ). 
CO fluxes are derived from the present observations, and H2 masses use a CO-to-H2 conversion factor of 3.0x10^" 
cm~^ (Kkms~^)~^. 



These observations were carried out in the C configura- 
tion (projected baselines 3 to 34 kA) in the spring and 
faU of 2003. Additional data for NGC 3032 and NGC 
4526 were obtained in the B configuration (projected 
baselines 4 to 83 kA) in March 2003. Total observing 
times were 8 to 17 hours per galaxy in the C configura- 
tion and 10 to 16 hours in the B configuration. A single 
pointing centered on the optical position was used for all 
galaxies. Each observation covered a velocity range of 
1300 km s~^, and these data have sensitivity to struc- 
tures from point sources up to objects 60-80" in diam- 
eter. System temperatures were mostly in the 300-500 
K range. Table [5] summarizes important parameters of 
these observations. 

Reduction of the BIMA data was carried out 

using standard tasks in the MIRIAD package 

(jSault. Teuben. fc Wright 1 19951 ). Electrical line length 
calibration was applied to all tracks. The C configura- 
tion data for NGC 3032 and 60% of the C configuration 
data for NGC 4526 were also explicitly corrected for 
amplitude decorrelation on longer baselines using data 
from an atmos pher ic phase monitor and the M IRIAD 
task uvdecor (iLavl [1999: Akeson 1998; Rega n et al.l 
I2OOII : IWongj |2001[ ). The atmospheric decorrelation is 
estimated using a small interferometer with a fixed 100 
meter baseline which measures the rms path length 
difference in the signal from a commercial broadcast 
satellite. Data requiring the decorrelation correction 
had rms path length differences in the range of 300 to 
700 microns, and the median amplitude correction factor 
was about 13%. However, the longest baselines in these 
datasets were later flagged because of their very poor 
phase stability. The remainder of the tracks were taken 
in stabler weather and were not explicitly corrected for 
decorrelation, because norr nal amplitud e calibration can 
take out most of the effect (|Wongll200lD . 

Absolute flux calibration was based on observations of 
Mars and 3C273, and comparisons of flux measurements 
on all the observed calibrators suggest that the absolute 
flux uncertainties are in the range of 15%-20%. Phase 
drifts as a function of time were corrected by means of a 
nearby calibrator observed every 30 to 40 minutes. Gain 



variations as a function of frequency were corrected by 
the online passband calibration system. Inspection of 
data for 3C273 indicate that residual passband variations 
are on the order of 10% or less in amplitude and 2° in 
phase across the entire band. 

The calibrated visibility data were weighted by the 
inverse square of the system temperature and the in- 
verse square of the amplitude decorrelation correction 
factor (if used), then Fourier transformed. Dirty images 
were lightly deconvolved with the Clark clean algorithm 
(|Clarklll980f ). as appropriate for these compact, rather 
low signal-to-noise detections. No continuum subtrac- 
tion was needed. Integrated intensity and mean velocity 
maps were produced by the masking method: the decon- 
volved image cube was smoothed along both spatial and 
velocity axes, and the smoothed cube was clipped at ±2.5 
times the rms noise in a channel. The clipped version of 
the smoothed cube was then used as a mask to define 
a three-dimensional volume in the original, unsmoothed 
cube in which the emission wa s integrated over velocity 
(|Wong|l2001l : lRegan et al.ll200l . Velocity maps were also 
constructed from Gaussian fits to the line profile at each 
position. 

Integrated spectra were constructed by first using the 
integrated intensity (moment 0) maps to define the spa- 
tial extent of the emission, then integrating over that 
fixed extent for every channel. Continuum images were 
also made by averaging all of the line-free channels in the 
final spectral line cubes, but no continuum emission was 
detected. Table [5] gives 3a limits for point source contin- 
uum emission at the centers of the galaxies, and Table [3] 
gives beam size and sensitivity information for the final 
spectral line cubes. 

4. RESULTS 

4.1. Total CO fluxes 

Comparisons of the total CO fluxes in the BIMA im- 
ages and in previous single dish data show that the in- 
terferometer usually recovers all of the CO emission. 
In the case of NGC 4459, the best single dish flux 
is from a survey of the SAURON early -type galaxies 
made with the IRAM 30m telescope by ICombes et al.l 
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TABLE 2 
CO Observation Parameters 



Galaxy 


Flux cal 


Phase cal 


Velocity Range 


FOV 


3mm cont. 








km 


kpc 


mJy 


NGC 3032 


3c273 


0927+390, 0854+201 


(876, 2128) 


10.5 


< 8 


NGC 4150 


3c273 


1159+292 


(-455, 766) 


6.6 


< 14 


NGC 4459 


3c273 


3c273 


(547, 1776) 


7.9 


< 11 


NGC 4526 


Mars, 3c273 


3c273 


(-233, 1050) 


7.9 


< 14 



Note. — Field of view (FOV) is the FWHM of the primary beam (100") at the distances in 
Table [1] 



TABLE 3 
CO Image Properties 



Galaxy 


Beam 


Beam 


Chan. 


noise 


N(H2) 




It 


kpc 


km s~^ 


Jy bm~^ 


10^0 cm-2 


NGC3032 


3.4x2.7 


0.4x0.3 


40.4 


0.013 


47.7 




8.1x6.7 


0.8x0.7 


20.2 


0.026 


8.0 




6.1x5.0 


0.6x0.5 


12.12 


0.025 


8.3 


NGC4526 


5.0x3.8 


0.4x0.3 


20.08 


0.021 


18.1 




4.3x3.1 


0.3x0.2 


40.16 


0.018 


45.5 


NGC4150 


8.5x5.1 


0.6x0.3 


40.04 


0.021 


15.9 




8.5x5.1 


0.6x0.3 


20.03 


0.029 


11.0 


NGC4459 


9.0x5.5 


0.7x0.4 


40.31 


0.014 


9.4 




9.0x5.5 


0.7x0.4 


20.14 


0.020 


6.7 



Note. — The N(H2) limit is a sensitivity estimate corresponding 
to a 3a- signal in one channel. 



(|2007l ). That work finds 54.0 Jykms-\ with a statis- 
tical uncertainty of 2.4 Jykms~^ and an absolute cali- 
bration uncertainty of 15% to 20%. The present BIMA 
data yield 56 zfc 11 Jykms~^. An ea rlier, much nois- 
ier detection by iThronson et~an ()1989D found 120 ± 55 
Jykms""'^. Thus, the single dish CO fluxes are en- 
tirely consistent with the flux measured from the in- 
terferometric map. We adopt here a CO-to-H2 conver- 
sion factor of 3.0x10^° c m~^ (Kkms~^)~^ and distances 
from iTonry et al.l (|2001h . scaled down by 3% following 
iMei et al.l (|2005i r At 15.7 Mpc, the H2 mass of NGC 
4459 is (1.6 ± 0.3) xlO^ Mq. 

Our BIMA-derived CO flux of NGC 3032 (93 ± 18 
Jykms^^) is also consiste nt with most of t he pr evious 
single dish measurements. IThronson et al.l ()1989D mea- 
sured 85 ± 25 Jyk ms"^ in the 45" beam of the FCRAO 
14m telescope, and ISage fc Wrobell (l989) measured 92 
± 10 Jykms-i in the 55" beam of the NRAO 12m 
Howev er, the more recent measurements of lCombes et al.l 
()2007l ) give only 46. 0± 1.5 Jykms"^ (again, with typi- 
cally 15% to 20% absolute calibration uncertainty) in the 
22" beam of the IRAM 30m telescope. Since the BIMA 
data show emission over a region at least 25" in diameter, 
it is plausible that some emission was missed by the 30m 
beam. At a distance of 21.4 Mpc, our flux corresponds 
to (5.0 ± 1.0) xlO^ Mq of H2. 

In contrast, the CO flux we flnd for NGC 4526 (180 
± 36 Jykms"^) is signi ficantly larger tha n previous sin- 
gle dish measurements. iLees et a Lldl SSLlD quote a value 
of 118±25 Jy kms-i measured bv'Sage fc Wrobell (|1989h 



with the NRAO 12m telescope, and lCombes et al.l ()2007f ) 
measured 118 ±4 Jykms"^ (statistical) at the IRAM 
30m telescope. Since the CO images show emission over 
a region nearly 30" in diameter, the 30m telescope may 
have missed significant flux. Furthermore, the spectrum 



of lSage fc Wrobell (|1989( ) is quite asymmetric, which may 
indicate mispointing or an inaccurate spectral baseline 
(neither of which plague the interferometric fluxes). For 
a distance of 16.4 Mpc, our flux corresponds to a mass 
of (5.7 ± 1.1) xlO* Mq of H2. Thus, in NGC 3032, 
NGC 4459, and NGC 4526 the evidence suggests that 
the interferometric maps have recovered all of the CO 
emission. 

There is some inconsistency in the single dish CO 
fluxes of NGC 4150, which makes it more difficult to 
judge whether the interferometric data are missing flux. 
The most recent 30m CO flux of NGC 4150 is 30.2 
zfc 2.4 (statis t ical) ib 5 (calibration) Jyknis"^ from 
ICombes etall ()2007( ). whereas the BIMA data yield 26 
± 5 Jykm s"^. These tw o mea surements are in good 
agreement. IWelch fc Sagi (|2003l ) measured 45± 2 ± 7 
Jykms"^, also with the IRAM 30m telescope, and 
this value may not be inconsistent with the previous 
tw o. However , a flu x of 77 ± 14 Jykms"^ reported 
bv lLerov et all ()2005l ) from the ARO 12m telescope (the 
old NRAO 12m telescope) is significantly higher, which 
might be due to an absolute calibration uncertainty, an 
underestimated baseline level, or even some extended 
molecular gas which is sampled by the 55" beam of the 
12m but not by the 22" beam of the 30m. At a distance 
of 13.4 Mpc, our BIMA flux corresponds to an H2 mass 
of (5.5 ± 1.1) xlO^ Mq. 

4.2. NGC 3032 

The CO in NGC 3032 is found in a centrally concen- 
trated structure of rather narrow linewidth. Emission 
is detected over 145 km s~^ centered on a systemic ve- 
locity of 1555 km s~^ (Figure [1]), so the CO systemic 
velocity is in good agreement with the stellar absorp- 
tion line velocity measu rements of 1555 ± 41 km s~^ 
from iFalco et al. (I1999D and 1559 ± 10 km s'^ from 
lEmsellem et alj ()2004[) . The diameter of the CO emis- 
sion is approximately 30", with a peak on the nucleus 
of the galaxy and a tail extending 15" to the southeast 
(Figure [2]). Based on a comparison with dust maps (see 
below), we infer that the molecular gas lies in a circular 
disk of radius 14" (1.5 kpc), yielding an average sur- 
face density of « 100 M© pc~^ including helium. NGC 
3032 is also known t o contain 0.9 Jykms"^ ( 1.0x10^ 
Mq) of HI emission (iDuprie fc Schneideil[l996l) . giving 
M(H2)/M(HI) « 5. 

Individual CO channel maps (Figure [3]) show good 
agreement between the distribution of molecular gas and 
dust in NGC 3032; dust and gas are both found in an 
inclined, rotating disk. An unsharp-masked WFPC2 im- 
age (Figure ^ shows a bright nucleus surrounded by a 
dark dusty ring (or two tightly wrapped arms) roughly 
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Fig. 1. — CO spectrum of NGC 3032. The spectrum was con- 
structed by first using the integrated intensity image (Figure [2J 
to define an irregular mask region within which the emission is lo- 
cated. The intensity was integrated over the same spatial region for 
every channel, so the noise in the line-free regions of the spectrum 
should be indicative of the noise on the line as well. 




RA (J2000) 



Fig. 2.— CO distribution in NGC 3032. The greyscale and fine 
black contours show the SDSS g image, with a contour interval of 
0.75 mag (a factor of two). White contours show the integrated CO 
distribution at —10, 10, 20, 30, 50, 70, and 90 percent of the peak 
(18.1 Jyb-ikms-l= 1.6x10^^ cm-^ or 260 Mq pc'^; helium is 
not included in the surface mass density). The hashed ellipse in 
the upper right corner shows the CO resolution (6.1" X 5.0"). 

6" in major axis diameter, and beyond that ring is a disk 
of flocculcnt spiral dust features interspersed with bright 
point sources. The dust features cover a region 28" in 
diameter. CO channels have a butterfly-wing structure 
with compact emission in the end channels and central 
channels elongated in the direction of the kinematic mi- 
nor axis. The extent of the CO emission matches that of 
the dust disk and the gas kinematic major axis matches 
the dust morphological major axis. The southeast "tail" 
of emission in Figure [5] is also visible in the channel maps 
at 1580 to 1616 km as emission at about the 75 
mjy beam~^ level, smoothly tracking to smaller radii 
with increasing velocity; it does not appear to be as- 
sociated with a dust feature. 



Kinematic analysis of the CO in NGC 3032 was made 
with several different techniques whose results are all in 
good agreement. To the velocity field shown in Figure [?] 
we fitted a model described by an inclined disk with the 
rotation curve V{r) oc 1 — e~^^''°, as encoded in the Na- 
tional Radio Astronomy Observatory's AIPS task GAL. 
The details of the shape of this rotation curve are not 
critical, but it has the desired behavior of rising quickly 
at small radii and asymptotically flattening. We also em - 
ployed the kinemetry analysis of Kra inovic et al.l (|2006) , 
and we made a tilted ring analysis with the rotcur task 
of the GIPSY package from the Kapteyn Institute, Ri- 
jksuniversiteit Groningen. In usi ng the rotcur task w e 
followed the procedure outlined bv lSwaters et al.l (|1999l ). 

The fitted kinematic center positions are within an arc- 
second of the optical nucleus measured from the SDSS 
g image. The systemic velocity of the galaxy is also ro- 
bustly fitted to be 1555 ± 1 km s~^. The inclination of 
the disk is poorly constrained by the CO data, however, 
and is more accurately taken from the axis ratio of the 
dust disk which gives 44° ± 4°. The Jeans dynamical 
modeling we did in this work is consistent and suggests 
i = 42°. 

Measurements of the globally averaged CO kinematic 
position angle fall in the range 90° to 97°, measured to 
the receding major axis, while the tilted ring model and 
the kinemetry both suggest a gentle twist from w 81° in 
the central resolution element to 97° at 8"< r < 14". 
This twist is also suggested by the hint of an integral- 
sign shape in the zero velocity curve of the velocity field, 
although warping is not obvious in the images of the 
dust disk. A major axis position-velocity diagram sliced 
through the data cube at 95° is shown in Figure [5] This 
position- velocity diagram shows the typical pattern with 
a steeply rising inner portion followed by a flattening 
(especially noticeable on the eastern side of the galaxy) 
beyond radii « 4". 

The kinematic position angle for the molecular disk is 
consistent with the photometric position angle, defined 
as that of the ellipse of inertia of the surface bright- 
ness (99.6°, CappeUari et al. 2007), bu t is 180° offset 
from the stellar kinematic an gle of —89° (|Emsellem et alj 
12004 ICappellari et all [20071 ).^ Thus, the molecular gas 
in NGC 3032 is counterrotating with re spect to the bulk 
of the stars in the galaxy. However, iMcDermid et al.l 
()2006airbh found a large radial age gradient in the galaxy 
and a counterrotating stellar core, which suggests recent 
star formation in the molecular disk. The sense of rota- 
tion of the molecular gas is consistent with that of the 
young counterrotating stars but inconsistent with that of 
the older, more extended population. 

4.3. NGC 4150 

NGC 4150 shows a relatively small amount of molec- 
ular gas in a very compact structure at the center of 
the galaxy. We fit the integrated spectrum in Figure 
[6] with a Gaussian whose central velocity is 239 ± 20 
km s~^, which we take as the systemic velocity of the CO 
emission. This CO velocity is in good agreement with 

^ ICappellari et aLl II2007I) quote a stellar kinematic angle of 
PAkin = 91° I they define it as it as the direction along which IV] 
is maximum. The velocity fields shown in EmscUem ct al. (2003) 
make it clear that this is to the stellar major axis on the approach- 
ing side rather than the receding side. 
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Fig. 3. — CO channel maps for NGC 3032. The velocity of each channel is indicated in the top left corner and the beam size in the 
bottom right corner of the first panel. Contour levels are —3, —2, 2, 3, 5, 7, 9, and 11 times the rms noise (25 mjy beam~^). The greyscale 
is an unsharp-masked WFPC2 image in the F606W filter. 



the stellar v elocity measurements of 226 ± 22 km s ^ 
(iFisher et al.llT995fl . 20 8 ± 30 km s'^ (iFalco et alJ ll999h. 
and 219 ± 10 km s'^ (jEmsellem et al.ll2004[ ). CO emis- 
sion is detected over a velocity range of 180 km s~^, and 
the bulk of this gas is within a few arcseconds of the 
nucleus of the galaxy (Figure [7]). 

Several dust structures are apparent in the center of 
NGC 4150 (Figure [8]). A dark dust lane bisects the nu- 
cleus from northwest to southeast, and the galaxy has 
more dust clouds a few arcseconds northwest and east of 



the nucleus. There is an irregular dust ring, stronger in 
the north, of semimajor and semiminor axes 7.5" x 4", 
with a fainter dust arm curving around the southeast side 
of the nucleus to radii « 12". The peaks in the channel 
maps suggest that most of the molecular gas is located in 
the bisecting dust lane, and the CO kinematic major axis 
roughly matches the major axis of the dust ring. Figure 
[7] shows a faint tail of emission stretching to 30" south 
of the nucleus, also visi ble in the channel map s from 265 
km s-i to 325 km s"^ iMorganti et all (|2006l) show that 
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Fig. 4.— CO velocity field of NGC 3032. The velocities are 
derived from Gaussian fits to the line profile at each position. Con- 
tours range from 1495 km s~^ to 1615 km s~^ at intervals of 20 
km s~^. The beam size is indicated in the top right corner. 




RA offset (arcsec; J2000) 

Fig. 5. — Major axis position-velocity diagram for NGC 3032. 
The galaxy is sliced through the center at a position angle of 95° . 
The abcissa is measured along the slice direction with the east side 
of the galaxy in positive values on the left side of the plot. The 
greyscale and black and white contours show CO intensity (con- 
tours are —3, —2, 2, 3, 5, 7, 9, and 11 times the rms noise). Red 
crosses show the mean stellar velocity from the SAURON data, 
green crosses show H/3 velocities and blue show [O III], The two 
solid magenta lines are two different estimates of the circular ve- 
locity Vcsini (see Section l5.3l l. 

the 2.5x10^ Mq of HI in the galaxy (M(H2)/M(HI) = 
23) is elongated roughly along the optical major axis in a 
structure of 1' diameter, with additional HI to the south 
and southwest, so the similarities between the CO and 
HI lend credence to the reality of the CO tail. 

Due to the small angular size of the CO emission the 
velocity field (Figure[9]) does not offer go od constraints on 
the ki nematic parameters of the galaxy. ICappellari et all 
()2007| ) quote both the photometric and stellar kinematic 
position angles as —33°, with uncertainties of a few de- 
grees. The orientation of the dust ring (Figure [H]) is 
—34° ± 5°, and although the CO velocity field is com- 
plex, it is not inconsistent with this value. A slice at po- 
sition angle —34° gives the major axis position-velocity 
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Fig. 6.— CO spectrum of NGC 4150. 
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Fig. 7. — CO distribution in NGC 4150. The greyscale and fine 
black contours show the SDSS g image, with a contour interval of 
0.75 mag (a factor of two). White contours show the integrated CO 
distribution at —10, 10, 20, 30, 50, 70, and 90 percent of the peak 
(17.0 Jyb-lkms-l= 1.1x10^^ ^^-2 helium is 

not included in the surface mass density). The hashed ellipse in 
the upper right corner shows the CO resolution (8.5" X 5.1"). 

diagram of Figure [TOl It appears that very little of the 
molecular gas is located beyond the turnover point in the 
rotation curve. The inclination of the galaxy is measured 
from the axis ratio of the dust ring to be 54°±5°, consis- 
tent with the valu e of 52° from the dynamical modelling 
(jCappellari et al.l[200& ). 

There is a counterrota ting stellar core in the central 
arcseconds of NGC 4150 (jMcDermid et al.ll2006allbl . and 
it is also visible in the position- velocity diagram of Figure 
[T0|) . This counterrotating core is believed to have formed 
recently as the galaxy shows a very strong radial age gra- 
dient. Curiously, though, the core does not appear to 
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Fig. 8. — CO channel maps for NGC 4150. The velocity of each channel is indicated in the top left corner and the beam size is in the 
bottom left corner of the first panel. Contour levels are —3, —2, 2, 3, 4, 5, and 6 times 29 mjybeam^^ (fa Icr). The greyscale is an F606W 
- F814W {oi R- I) image constructed from HST WFPC2 data. 



have a counterpart in the molecular gas at the current 
resolution and sensitivity. The counterrotating core is 
also not obvious in the [O III] velocities of NGC 4150 
(Figure [TO]) , and indeed the CO velocities of this galaxy 
are a better match to the [O III] velocities than they are 
to the stellar velocities. A good match between CO ve- 
locities and ionized gas velocities suggests the possibility 
that the ionized gas could trace star formation activity, 
but the relationship to the counterrotating stellar core is 
not yet clear. 



4.4. NGC 4459 

The CO emission from NGC 4459 is also in a central 
disk, with hints of a double-horned structure in the in- 
tegrated spectrum (Figure [TT]). A double-horned spec- 
trum is typically produced by gas in the flat part of 
the rotation curve. Emission is detected over a velocity 
range of 400 km s^^ centered near a systemic velocity 
of 1210 ± 20 km s^^ (but see below). The CO veloc- 
ity is thus in good agr eement with the st ellar velocity of 
1232 ± 40 km s^^ from lFalco et all ()1999f ) and 1200 ± 10 
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Fig. 9.— CO velocity field of NGC 4150. The velocities shown 
here are the intensity- weighted mean velocity (moment 1) rather 
than Gaussian fits, as the fits to individual profiles are not robust 
in this rather low signal-to-noise detection. Contours are from 180 
km s~^ to 320 km at intervals of 20 km s~^. The beam size is 
indicated in the top left corner. 
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Fig. 10. — Major axis position-velocity diagram for NGC 4150. 
The galaxy is sliced through the center at a position angle of —34°. 
The abcissa is measured along the slice with the southern side of the 
galaxy in negative values. The greyscale, black and white contours 
show CO intensity (contours are —2, 2, 3, and 5 times the rms 
noise level). Red crosses show the mean stellar velocities from the 
SAURON data and blue crosses show [O III]. (Velocities of H/3 
are not fitted independently of [O III] in this galaxy.) The solid 
magenta line is Vcsini. 

km s^^ from lEmsellem et al.l (|2004D . 

The integrated intensity map (Figure fT^ shows a com- 
pact structure centered on the galaxy nucleus, no more 
than 20" in diameter. Little structure is evident, and 
there is no strong evidence for arms or tails of molecular 
gas. The individual channel maps in Figure [13] again 
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Fig. 11.— CO spectrum of NGC 4459. 
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Fig. 12.— CO distribution in NGC 4459. The greyscale and fine 
black contours show the SDSS g image, with a contour interval of 
0.75 mag (a factor of two). White contours show the integrated 
CO distribution at -10, 10, 20, 30, 50, 70, and 90 percent of the 
peak (21.4 Jy b-^ kms-i= 1.2x10^2 cm-^ or 190 Mq pc-^). The 
hashed ellipse in the upper right corner shows the CO resolution 
(9.0" X 5.5"). 

show a typical disk pattern with higher intensities at 
the extreme velocities than near the systemic velocity. 
The centroids of the emission in each channel show ex- 
cellent agreement with the extent of the prominent, well- 
developed dust disk. The dust fills a flocculent disk with 
a rather sharp outer edge at a semimajor axis of 8.5" 
and only a few faint spiral dust features beyond. From a 
radius of about 3.75", two prominent dusty spiral arms 
can be traced inward to an inner ring of semimajor axis 
2.0". The larger disk has a major axis position angle of 
102° but the inner ring appears to have a slightly dif- 
ferent orientation, with position angle « 90°. The bulk 
of the CO emission appears associated with the larger- 
scale, 8.5" disk. For a circular disk of radius 8.5" (670 
pc), the average gas surface density is 170 IVIq pc^^, in- 
cluding helium. The north side of the disk must be the 
near side, suggesting that both the flocculent and inner 
spiral arms are trailing. 
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Fig. 13. — CO channel maps for NGC 4459. The velocity of each channel is indicated in the top left corner and the beam size in the 
bottom right corner of the first panel. Contour levels are —3, —2, 2, 3, 5, 7, and 9 times 14 mjybeam^^ (fs Icr). The greyscale is the 
inner portion of an unsharp masked F555W image from WFPC2. 

analysis of the stellar velocities gives position angles of 
280° outside the dust disk, decreasing to 274° inside it 
(|Krainovic et al.ll2006r ). A global kinematic position an- 
gl e measurement for the m olecular gas, measured both as 
in lKrainovic et all (|2006l ) and by fitting a model velocity 
field, yields values in the range 269° to 273° to the reced- 
ing major axis. The molecular gas kinematics are thus 
very well aligned with the stellar kinematics and the stel- 
lar isophotes. A tilted ring analysis suggests the systemic 
velocity to be 1196±2 km s'^. The major axis position- 



The velocity field (Figure [T4|) shows predominantly 
straight, parallel isovelocity contours, so again the in- 
clination of the gas disk is better obtained from the dust 
disk than from the gas itself. The axial ratio of the dust 
disk gives an estimated inclination of 46 ° ±2°, consistent 
with 47° from dynamical modelling (jCappe llari ct al. 
l2006f ). For the kinematic analysis we also constrain the 
kinematic center of the galaxy to be the position of the 
optical nucleus. Stellar isophotes have a position an- 
gle of 280° outside the dust disk, while a kinemetric 
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Fig. 14.— CO velocity field of NGC 4459. Velocities are derived 
from Gaussian fits to the line profile at each position. Contours 
are from 1050 km s~^ to 1350 km s~^ at intervals of 50 km s~^. 
The beam size is indicated in the top left corner. 
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Fig. 15. — Major axis position-velocity diagram for NGC 4459. 
The galaxy is sliced through the center at a position angle of 90° . 
The east side of the galaxy is at the left side of the figure. The 
greyscale, black, and white contours show CO intensity (contours 
are —3, —2, 2, 3, 5, 7, and 9 times the rms noise level). Red 
symbols show the SAURON stellar velocities along the major axis; 
blue symbols show [O III] velocities and green symbols are H/3. 
The solid magenta line is Vc sini. 

velocity diagram (Figure llSp is very strongly peaked at 
its maximum velocities, suggesting that the bulk of the 
molecular gas is located in the flat part of the galaxy's 
rotation curve. 

4.5. NGC 4526 

In NGC 4526 the CO lies in a well-developed, nearly 
edge-on disk. The integrated CO profile shows the famil- 
iar double- horned shape fFigure fTB)) . with both horns ap- 
proximately equally bright. The systemic velocity of the 
galaxy is 613 ± 10 km s^^ and emission is detected over 
a total width of 683 km s^^. A line width of this mag- 
nitude, while not unknown, is unusually large for lentic- 



ulars. In fact, we argue in section [5731 that the galaxy's 
maximum circular velocity is about 355 km s"^ and in 
the compilation of 243 galaxies (54 lenticulars) made by 
ICourteau et al.l (|2007| ). NGC 4526's circular velocity is 
exceeded by only two other lenticulars, one Sa galaxy, 
and 10 ellipticals. 

The observed CO systemic velocity is consistent with 
the 



stellar velocit y 
Falco et al.l |1999D_ 



of 592 ± 48 km s~^ given by 
an d 626 ± 10 km s~^ given by 
Emsellem et al.l ( 2004 ). but is inconsistent with the HI 
velocity of 448 ± 8 km s~^ (jde Vaucouleurs et al]ll991t 
iDavies fc Lewiill973[ ). In fact, a com parison of Figure 
[T6]with the HI spectrum in lDavies fc L ewis (1973) shows 
that the latter authors only detected one "horn" , so their 
HI velocity is not a good estimate of the systemic veloc- 
ity. Interferometric HI observations would undoubtedly 
be useful in revealing the relationships between atomic 
and molecular gas. 

The CO distribution (Figure [T7|) is elongated in the 
direction of the optical major axis. Fitting a two- 
dimensional Gaussian to the integrated intensity image 
reveals that the disk is poorly resolved in the short di- 
mension; the deconvolved minor axis FWHM is 4.3", 
comparable to the beam size. The position angle of this 
gas disk is estimated as —72.7° ± 1.8°. 

Individual channel maps (Figure [TH]) show a close cor- 
respondence between the CO and dust disks. Unsharp- 
masked HST images show that the dust is confined to 
radii less than « 14" - 15". Although the inclination 
is quite high, there appear to be at least two dominant 
dust rings of radii « 14" and 10". (These radii are 
uncertain because the features are most prominent along 
the minor axis.) The centroids of CO emission in the ex- 
treme channels suggest that the bulk of the CO may be 
associated with the inner of these two rings. Assuming 
the gas to be in a circular disk of radius 14" (1.1 kpc), the 
average surface density is 200 Mq pc~^ including helium. 

As in the case of NGC 4459, the velocity field (Figure 
[T9|) shows straight, parallel isovelocity contours through- 
out; it constrains the kinematic position angle but not 
the inclination of the gas disk. However, the dynamical 
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Fig. 17. — CO distribution in NGC 4526. The greyscale and fine 
black contours show the SDSS g image with a contour interval of 
0.75 mag (a factor of two). White contours show the integrated 
CO distribution at -10, 10, 20, 30, 50, 70, and 90 percent of the 
peak (37.1 Jyb-lkms-l= 5.3x10^2 cm'^ or 850 Mq pc'^). The 
hashed ellipse in the upper right corner shows the CO resolution 
(5.0" X 3.9"). A slightly higher resolution integrated intensity im- 
age has a peak of 32.2 Jyb-^kms"! = 6.8x10^^ cm-^ or 1080 
Mq pc~2 at a beam size of 4.3" X 3.1". 

modelling of lCappellari et all (|2006l ) gives an inclination 
of 79° which is consistent with values of 75° ± 2° inferred 
from the axis ratio of the dust disk. Fits of a model ex- 
ponential velocity field show that the kinematic center of 
the molecular gas is consistent with the position of the 
optical nucleus, within 0.5". The global CO kinematic 
position angle, calculated from those fits and also via 
the method of lKrainovic et al.l (|2006f) . is -78° ± 3° to 
the receding major axis. More detailed kinemetric anal- 
ysis suggests a gradual trend from PA —78° in the inner 
resolution element to —72° at radii > 6". In compari- 
son, the stellar photometric major axis is —67° and the 
stella r kinematic major axis is -69° (iCappell ari et al.l 
l2007f) , so the molecular gas (especially at its outer radii) 
is aligned with the stellar body of the galaxy to within a 
few degrees. Figure 1201 the major axis position- velocity 
diagram, clearly shows that the gas traces a turnover in 
the rotation curve with concentrations of gas near the 
turnover radius. There is a modest degree of asymme- 
try, with somewhat stronger CO emission on the receding 
(west) side of the major axis. 

5. DISCUSSION 
5.1. Origins of the Molecular Gas 

In NGC 4459 and NGC 4526 the molecular gas kine- 
matics are consistent with both the stellar photometric 
axes and the stellar kinematic axes. The differences are 
on the order of a few degrees, which is the level of pre- 
cision available from existing data. In NGC 4150 the 
molecular gas kinematics are more difficult to describe 
but the orientation of the dust ring is also consistent 
with the stellar photometric and kinematic axes to a de- 
gree or better. In NGC 3032 the molecular gas kinematic 
position angle is consistent with the photometric axis to 
within a few degrees, but it is nearly 180° away from the 
stellar kinematic axis. Despite the dramatic counterrota- 
tion of the gas in NGC 3032, these alignments imply that 
all four galaxies have nearly axisymmetric potentials and 
that the gas is well settled into the equatorial plane. 

In NGC 4150, NGC 4459, and NGC 4526 the agree- 
ments between the sense of rotation of the molecular gas 



and stars also imply that the molecular gas could have 
originated in internal stellar mass loss. Such an internal 
origin is not required by the data, of course, as the gas 
could have been captured into prograde rotation from 
outside sources, and internal secular evolution over sev- 
eral dynamical timescales would gradually bring it to the 
galaxy's equatorial plane. (Orbital periods are a few 10^ 
yr at the edges of the disks.) 

In contrast to the other three galaxies, the conspic- 
uous counterrotation of the molecular gas and stars in 
NGC 3032 provides clear evidence that this molecular 
gas could not have originated in internal stellar mass 
loss. The only scenario that could explain the coun- 
terrotation of internally produced gas is to invoke some 
gravitational interaction which strongly torques the gas 
(possibly at large radii) while leaving the remainder of 
the galaxy seemingly untouched. This scenario seems un- 
necessarily complex, however, so we suggest instead that 
the molecular gas was captured through "cold" accretion 
from the intergalactic medium or in an interaction or a 
minor merger with a gas-rich neighbor. Perhaps it is even 
a remnant of a major merger which formed the present 
galaxy. The high degree of regularity in the gas kinemat- 
ics and stellar morphology suggests that this event did 
not occur recently, as the orbital period at the outer edge 
of the CO disk is on the order of 10^ yr. Furthermore, 
modeling the optical image (Section 15. 2|) shows that the 
isophotes are very regular to at least 70" (7.3 kpc, or five 
times the radius of the molecular gas). 

NGC 3032 is significant, then, because it is one of the 
few early-type galaxies in which the possibility of an in- 
ternal origin for the molecular gas can be fi rmly excluded. 
It con tradicts the picture outlined by ISage fc WelchI 
()2006l ). who suggested that most of the molecular gas 
in lenticulars should have an internal origin like that de- 
scribed by Temi, Brightenti, fc Mathews (2007). The ev- 
idence for the suggestion of ISage fc WelchI (|2006l ) came 
from contrasting the atomic and molecular properties of 
lenticulars: the atomic gas, which is usually more ex- 
tended and shows different kinematics than the molec- 
ular gas, tends to dominate the gas content in gas-rich 
galaxies, whereas th e molecular phase dom inates in gas- 
poor galaxies. Thus ISage fc WelchI (j2006f ) hypothesized 
that the atomic gas could be primarily attributed to ex- 
ternal sources but the molecular gas to internal sources. 
While that picture may still be roughly correct, it cannot 
account for the molecular gas in NGC 3032. 

On the other hand, it must still be true that the evolved 
stars in NGC 3032 exp elled mass into the ISM . Ac- 
cording to the mo dels of iFaber fc Gallagheil ()1976f ) and 
iCiotti et al.l(fT99lh . given the galaxy's present luminosity 
the stellar mass loss would amount to 4x10* to 6x10^ 
Mq over a Hubble time. We have not detected that gas 
yet, as it would show prograde rotation. The stellar 
mass loss must now be in the form of hot gas, or per- 
haps it has been removed from the galaxy entirely (e.g. 
iTemi. Brightenti. fc Mathewsl[2007l ). 

It is not unusual to find counterrotating ionized gas 
in lenticular galaxies. Fro m their own data and a com- 
pilation of the hterature, iBureau fc Chund ()2006[ ) find 
that 15% ±4% of all SOs contain counterrotating ion- 
ized gas while 23% ±5% of the SOs with ionized gas 
have that ga s in re trogra de motion. Earlier work by 
iBertola et all (|1992l ) and iKuijken. Fisher, fc Merrifieldl 
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Fig. 18. — CO channel maps for NGC 4526. The velocity of each channel is indicated in the top left corner and the beam size in the 
bottom right corner of the first panel. Contour levels are —3, —2, 2, 3, 5, 7, and 9 times 19 mjy beam~^ (sb la). The greyscale is an 
unsharp-masked F555W image from WFPC2. 



is in aereement with this result. A more de- 
tailed study o f stellar - ionized gas misalignments was 
carried out by ISarzi" "eTall ([2005) using integral-field ob- 
servations. They find that among the lenticulars of the 
SAURON early-type sample, 9 of 20 galaxies have a star- 
gas kinematic misalignment greater than 30°. NGC 3032 
is the only one of the 9 that has a misalignment > 150°, 
though, so it is somewhat unusual in having ionized gas 
so nearly retrograde. We should be careful not to make 
assumptions about the behavior of the molecular gas 
based on the behavior of the ionized gas, however, since 
in general it is not obvious that these two phases of the 
interstellar medium are closely linked. For example, in 



NGC 4150 the ionized gas emission (ISarzi "eiTaH [20051 ) is 
distributed over a much larger scale than the molecular 
gas; ionized gas fills the 30" x 40" SAURON field of view 
whereas molecular gas is concentrated in the central 6". 
Clearly, a larger sample of resolved molecular maps of 
lenticular galaxies will be required before we can make 
better judgments about the origin of their molecular gas. 

5.2. Deriving circular velocities 

Molecular gas is naturally cold and dissipational and 
therefore, in equilibrium in an axisymmetric potential, 
it should settle in circular orbits in the equatorial plane. 
To the extent that the velocity dispersion of the molecu- 
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Fig. 19.— CO velocity field of NGC 4526. The velocities are 
derived from Gaussian fits to the line profile at each position. Con- 
tours are from 310 km s~^ to 910 km s~^ at intervals of 50 km s~^. 
The beam size is indicated in the top left corner. 
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Fig. 20. — Major axis position-velocity diagram for NGC 4526. 
The galaxy is sliced through the center at a position angle of 111°. 
The east side of the galaxy is at the left side of the figure. The 
greyscale, black, and white contours show CO intensity (contours 
are —2, 2, 3, 5, 7, and 9 times the rms noise level). Red symbols 
show the SAURON stellar velocities along the major axis; blue 
symbols show [O III] velocities and green symbols are H/3. The 
solid magenta line is Yc sini. 

lar gas is low (typically on the order of 10 to 20 km s^^) 
it should then provide an excellent tracer of the circu- 
lar velocity of the galaxy. And the circular velocity is, 
ultimately, the dynamical indicator of the total matter 
content of a galaxy. 

It is also possible to estimate the circular velocity 
from stellar kinematic data, with additional complica- 
tions. For a given edge-on axisymmctric potential, di- 
mensional arguments suggest that it may be possible 
to uniquely recover the three-dimensional (3D) orbital 
distribution of the stars in a galaxy with another 3D 
quantity, namely the knowledge of the stellar line-of-sight 
velocity-distribution at every position on the galaxy, as 
can be obtained with integral-field spectroscopy. It seems 
however unlikely that the axisymmctric potential itself, 
which is another two-dimensional function can also be 
uniquely re covered from the sa me 3D observations (see 
Section 3 of iValluri et al.|[200^ . for a discussion). In the 
stellar dynamical models either a constant mass-to-light 
ratio {M / L) or a parametric form for the potential are 
assumed. In this context it is worthwhile to compare the 



two different estimates of circular velocity in early-type 
g alaxies coming from diff erent assumptions. 

ICappellari et al.l (I2006D have us e d the stellar kinemat- 
ics presented in lEmsellem et al.l ()2004l ) to derive high 
quality dynamical M/L for a subset of the SAURON 
\ early-type galaxies, under the assumption of a constant 
3 M/L. High resolution and wide-field images of the galax- 
ies are modeled with a Multi-G aussian Expansion (MGE; 
lEmsellem et al.lll994t ICappe llari 2002) so that, if the in- 
clination of the galaxy (assumed to be axisymmctric) is 
known, the fitted MGE model of the surface density can 
be deprojected and a three-dimensional stellar distribu- 
tion can be estimated. Full two-dimensional maps of the 
stellar mean velocities, dispersions, and Gauss-Hermite 
hs-he parameters within about an effective radius are 
then used in combination with the stellar distribution 
to construct self-consistent two-integral Jeans and three- 
integral Schwarzschild models. The data are consistent 
with the assumption of a constant dynamical M/L within 
an effective radius, and in galaxies with regular dust disks 
the inclination inferred from the Jeans models is consis- 
tent with the axial ratio of the dust disk. 
Here we compute the circular velocities: 

where the gravitational potential ^(R ,z) is computed 
from the MGE models tabulated in ICappellari et al.l 
(|2006[ ). deprojected at their best fitting inclination, and 
scaled by the constant best-fitting Schwarzschild's M/L. 
For NGC 3032 no previous dynamical model exist so the 
MGE model, inclination and M/L were determined as 
described below. The circular velocity curves of the best- 
fit gravitational potentials for NGC 3032, NGC 4150, 
NGC 4459 and NGC 4526 are shown in Figures El [TOl 
[m and [201 

For NGC 3032 we constructed an MGE mo del for the 
photometry using the software of ICappellari (|200 2). To 
reduce the effect of dust extinction in the derivation of 
the stellar density distribution, in the MGE fit we com- 
bined an archival near-infrared HST/NICMOS/F160W 
(i7-band) image of the dusty central regions with a 
larger-field image taken with HST/WFPC2/F606W {V- 
band). The y-band image was adopted as photomet- 
ric reference and used to converted the MG E model 
to the Johnson band using the calibration of iDolphinI 
()2000f ). The distance-independent parameters of the 
PSF-deconvolved MGE models, c orrec ted for galactic ex- 
tinction following lSchlegel et al.l ()1998D . and adopting an 
absolute mag nitude of the Sun M (n.\/ = 4.83 mag (Ta- 
ble 2.1 of Binn ev fc Merrifieldlll998[ ).'are given in TableH 
We constructed a self-consistent two-integral axisym- 
mctric MGE Jeans model for the second velocity mo- 
ments /i2 = \nP + a^ of NGC 3032 and we determined 
the best-fitting (M/L)v ~ 1-7 and inclination i « 42°, as 
in Cappellari et al. (2006). Using the V - I = 1.07 color 
of iTonrv et all (|2001 ) this translates into {M/L)i ^ 1.3, 
which is among t he lowest measured v alues for any 
SAURON galaxy (jCappellari et al.ll2006[ ). This is stiU 
consistent within the scatter with their {M/L) — a re- 
lationship, consideri ng the low value of tX e = 90 km 
of NGC 3032 from EmscUra^eEaD ([200l . The low dy- 
namical M/L is also expected from their relation be- 
tween {M/L) - H/3 and the fact that NGC 3032 has the 
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TABLE 4 

MGE PARAMETERS FOR THE 
DECONVOLVED V-BAND STELLAR 
SURFACE BRIGHTNESS OF NGC 3032 



j 


log/; 


log <Tj 








(arcsec) 




1 


6.055 


-1.426 


0.895 


2 


5.117 


-0.775 


0.895 


3 


4.236 


-0.485 


0.895 


4 


3.882 


-0.019 


0.840 


5 


3.005 


0.282 


0.950 


6 


2.797 


0.597 


0.863 


7 


1.769 


0.913 


0.950 


8 


2.060 


1.098 


0.800 


9 


1.635 


1.532 


0.806 



Note. — is the peak surface 
brightness of the j-th Gaussian, which 
has dispersion aj and observed axial 
ratio q'j. 

lowest average H/3 line-strength among all the galaxies in 
the S AURON sample of E/SO galaxies (|Kuntschner et aLl 
l2006f ). This confirms that the variations in the stellar 
populations are the main driver for the observed differ- 
ences in the dynamical M/L. However the self-consistent 
model does not provide a good fit to the data, so the fit- 
ted M/L represents an average value for the central re- 
gions of the galaxy which are sampled by the SAURON 
kinematics. In particular the self-consistent two-integral 
model shows a much steeper radial decrease of ^2 along 
the major axis than observed. This is the opposite of 
what is generally measured in fiattened galaxies and disks 
(e.g.|Cappcllari et al. 2007.) and suggests that anisotropy 
cannot explain the discrepancy. An improved non-self- 
consistent model will be discussed in Section [5T3l 

5.3. Comparing circular velocities with the observed CO 

velocities 

We make a comparison between the stellar and the 
molecular indicators of circular velocity with tilted ring 
models, projected and convolved to the same resolution 
and sampling as the CO data. The GALMOD routine in 
the GIPSY software package is used to generate the mod- 
els. We specify the rotation velocities of the rings as com- 
puted in Section [521 and listed in Tables [3 [U and [71 The 
inclinations of the gas disks are also known, as described 
above, and the velocity dispersion in the molecular gas is 
assumed to be 10 km s~^. The exact value of the velocity 
dispersion used is not critical, as long as it is smaller than 
the channel width and consistent with the sharp cutoffs 
in emission at the outlying velocities. Nominally the ro- 
tation velocities should also be corrected for this finite 
velocity dispersion (the asymmetric drift effect), but as 
the rotation velocities are 100 to 350 km s~^ the correc- 
tion is negligible. The optical dust images also suggest a 
very thin molecular disk, here assigned a scale height of 
1" (again the exact value is not critical as it is smaller 
than the spatial resolution) . A gradual linear decrease in 
the gas surface density is assumed, with a sharp cutoff 
at the edge of the dust disk. The rotating disk model is 
then projected onto the plane of the sky at the specified 
inclination and the model emission is sampled at a pixel 
size and channel width matching the CO observations. 
The sampled cube is also spatially smoothed to have the 
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Fig. 21. — Observed and modeled position-velocity diagram for 
NGC 4526. The greysc ale, black, and white contours show the 
CO intensity as in Fig. 1201 red contours show the correspond- 
ing position-velocity diagram for the tilted ring model in which 
the ci rcular velo c ity is everywhere 5% smaller t han that given by 
Capp ellari et al] 1120071 ). as discussed in Section 15.31 The adjusted 
Vc sinj are shown in the blue line. Contour levels for the model 
position-velocity diagram are at the same levels (relative to the 
peak intensity) as for the data. 

same resolution as the CO beam. Finally the smoothed 
model cube is sliced along the galaxy major axis just as 
for the real data, and the slice is displayed with analo- 
gous contour levels (as percentages of the peak intensity, 
which is arbitrary). 

Figure [nH compares the data and the model position- 
velocity diagrams for NGC 4526. The assumed gas sm- 
face densities drop linearly by a factor of 5 between 1" 
and 14" and are zero beyond that. Placing the gas in 
rings rotating at the circular velocity produced a model 
with a bit too large velocity amplitude; the peaks in the 
model position-velocity diagram, tracing gas along the 
line of nodes, were separated by 703 km s~^ whereas in 
the data they are separated by 663 km s~^ (the difference 
is two channels). Therefore, Figure [21] actually shows a 
model in which the rotation velocity of each ring is 5% 
smaller than the circular velocities derived with the Al / L 
of iCappellari et al.l (|2006l ). and it shows that this model 
produces a very close match to both the ridgeline and the 
envelope of the observed emission. Even the rounded 
shoulder of the circular velocity curve (specifically, the 
relatively shallow but still rising slope in 1" < r < 8") 
is confirmed by the CO data. We have not made an 
exhaustive search in the rotation velocity or gas distri- 
bution parameter space, so we certainly cannot claim 
to have made the unique best match to the observed 
position-velocity diagram. However, the good agreement 
between the model and the data should be regarded as 
independent confirmat ion that the circu l ar vel ocities de- 
rived with the M/L of lCappellari et al.l (|2006l ) are accu- 
rate to about 5%. The M/L is thus accurate to about 
10%, roughly consistent with the quoted uncertainty. 

Figure [22I presents a similar comparison for NGC 4459. 
In this case the circu l ar vel ocities obtained with the M / L 
of ICappellari et all (|2006f ) are used without correction 



16 



TABLE 5 
Tilted Ring Models for NGC 
4526 









^gas 


It 


km s^^ 


km s^-"- 




1.34 


289.0 


274.6 


1.00 


1.80 


295.2 


280.4 


0.97 


2.42 


302.7 


287.6 


0.93 


3.26 


318.0 


302.1 


0.88 


4.37 


343.3 


326.1 


0.81 


5.88 


366.4 


348.1 


0.72 


7.90 


373.8 


355.1 


0.59 


10.62 


366.6 


348.3 


0.43 


14.26 


356.8 


339.0 


0.20 



Note. — The velocities Vrot(i) a-re 
equal to the stellar-derived circular 
velocities, and Vrot(2) = 0-95Vrot(i)- 



TABLE 6 
Tilted Ring Model for 
NGC 4459 



Radius 




Spas 


11 


km s"'^ 


1.5 


272.3 


1.00 


3.0 


268.7 


0.86 


4.5 


269.6 


0.73 


6.0 


266.8 


0.60 


7.5 


263.3 


0.47 


9.0 


262.1 


0.10 



and an excellent match is achieved, indicating again the 
high degree of accuracy of the dynamical M / L. Figures 
[T5]and[^also show that in NGC 4459 and NGC 4526 the 
molecular gas rotation speeds are 1.6 to 2.0 times larger 
than the mean stellar velocities at the same radii. In 
these galaxies the stars have significant pressure support, 
but the molecular kinematics and the stellar kinematics 
are both consistent with the same gravitational potential. 

A broader implication of this agreement is that 
the molecular and stellar kinematics for galaxies like 
NGC 4459 and NGC 4526 can be used with confidence in 
studies of th e circular velocity and the bulge velocity dis- 
persion fe.g. lCourteau et ani2007f) . where it is necessary 
to probe the relationships between the two dynamical 
indicators. For ex ample, several authors (e.g. iHol [20071 : 
[Shields et al.ll2006h have studied the development of the 
relationship between a galaxy's black hole mass and its 
bulge velocity dispersion. They have advocated use of 
the CO line width to trace the galaxy's maximum cir- 
cular velocity, hence its velocity dispersion (especially at 
high redshift where it is currently easier to measure the 
CO line width than other dynamical indicators of the 
galaxy mass) . These studies use a CO Tully-Fisher rela- 
tion to demonstrate that the CO line widths do indeed 
trace their host galaxies' maximum circular velocity. Our 
detailed comparison of molecular kinematics and circu- 
lar velocities supports this claim when the CO velocities 
and circular velocities at large radii are considered, even 
for early-type galaxies where the Tully-Fisher relation is 
not traditionally applied. 

Figure [TUl compares the CO position- velocity diagram 
in NGC 4150 to the inferred circular velocity. In this case 
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Fig. 22. — Observed and modeled position-velocity diagram for 
NGC 4459. The greyscale, black, and white contours show the 
CO intensity as in Fig. 1151 red contours show the model position- 
velocity diagram i n which the rotation vel ocity of each ring is the 
circular velocity of ICappellari et al.l I I2006I) (blue line). 

the total velocity range covered by the molecular gas is 
nearly as large as what would be expected from the cir- 
cular velocity curve. The CO velocities are also nearly 
equal to the mean stellar velocities at radii > 10", which 
is qualitatively consistent with the fact that the ob- 
served stellar velocity disper sions are much smaller than 
in NGC 4459 and NGC 4526 (jEmsellem et al.|[200l . It is 
worth noting that the adopted inclination for the galaxy, 
while consistent with the axis ratio of the dust ring, may 
not be appropriate for the molecular gas which is con- 
centrated in the nuclear dust lane. Detailed models are 
probably not yet useful for these data, and the most that 
can be said is that the inferred circular velocity curve is 
not inconsistent with the observed CO and stellar veloc- 
ities. 

In contrast to the cases of NGC 4459 and NGC 4526, 
there are significant disagreements between the inferred 
circular velocity and the observed molecular gas rotation 
speeds of NGC 3032. Assuming a constant dynamical 
M / L within an effective radius, the circular velocity rises 
quickly to a peak of 330 km s~^ at 0.08" and makes a 
sharp, nearly Keplerian decline through the inner arcsec- 
onds. It drops by a factor of two to 157 km s~^ at 5" 
and thereafter drops more slowly, passing through 132 
km s~^ at 10" which is nearly the edge of the molecular 
disk. But as Figure [5] shows, these velocities are much 
larger than the observed CO velocities especially interior 
to 5". 

The strong radial gradient i n the H/3 absorption lin e 
strength index of NGC 3032 (iKuntschner et al.ll2006f) 
and t he young inferred central age ([McDermid et al.l 
l2006a^ suggest that the stellar population might be 
better described with a variable M/L. Indeed, if the 
line s trength indice s are interp reted with the mod- 
els of iThomas et all (|2003[ ) and iMarastonI ()2005( ). the 
SAURON data are consistent with a decline of a fac- 
tor of two to three in the local stellar M/L from 10" 
to 1", even in the H band (from which the MGE models 
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Fig. 23. — Observed and modeled position- velocity diagram for 
NGC 3032. The greyscale, black, and white contours show the 
CO intensity as in Fig. [S] red contours show a model position- 
velocity diagram in which the rotation velocity of each ring is the 
circular velocity (blue line), assuming a decrease of a factor of two 
in the M/L interior to 10". The green line shows a rotation curve 
Vrot sini which reproduces the observed gas velocities, as in Table 

El 

are derived) . Moreover the inability of the self-consistent 
Jeans model of Section [5?2l to qualitatively reproduce the 
observed kinematics also indicates the need for a varying 
dynamical (total) M/L. Thus a second MGE model was 
constructed, with a variable M/L, such that the depro- 
jected Gaussian components having widths (dispersions) 
< 1" have M/L values half as large as the components 
with dispersions > 10", and the variation is logarith- 
mic between those extremes. Using this mass model a 
non-self-consistent Jeans model was computed and again 
fitted to the observed SAURON stellar kinematics to de- 
rive a new M/L. Interestingly this revised Jeans model, 
which has an M/L gradient dictated by the stellar pop- 
ulation measurements, now also well reproduces the ob- 
served SAURON stellar kinematics. The inferred circular 
velocity curve has a somewhat smaller peak velocity, 285 
km s^^ at 0.08", a steeper decline in the inner few arc- 
seconds and a shallower slope at 10", but its amplitude is 
the same at 10" (132 km s~^) as for the model with the 
constant M/L. Figure [23] shows the predicted CO ve- 
locities that would have been observed if the molecular 
gas in NGC 3032 were rotating at the inferred circular 
velocity (Table [7]). Clearly the assumption of a variable 
M/L does not remove the discrepancy between observed 
and expected velocities. 

A reasonably accurate measure of the discrepancy be- 
tween observed and expected velocities can be made from 
a model which describes the CO velocities well. Ta- 
ble [7] and Figure [53] show the rotation velocities in this 
model; they have been altered by hand to reproduce the 
observed position-velocity diagram. As before, no for- 
mal fitting has been done so the model shown cannot be 
claimed unique or best and its gross properties should be 
trusted to a greater degree than its details. This "work- 
ing" model has a rotation velocity of 98 km s""'^ at 6.5" 
and 113 km s^^ at 13.5". At the outer edge of the CO 
disk its velocity is only 10% smaller than that of the in- 



TABLE 7 

Tilted Ring Models for, NGC 
3032 



— — 

Radius 


— — 

' circ 






II 


km s~^ 


km s ^ 


1.5 


183.6 


88.7 


1.0 


2.7 


166.1 


89.8 


0.9 


3.6 


155.3 


90.9 


0.8 


4.5 


148.8 


92.0 


0.7 


5.5 


144.5 


93.2 


0.6 


6.5 


141.4 


97.6 


0.5 


7.5 


138.7 


101.0 


0.4 


8.7 


135.7 


106.3 


0.3 


10.0 


132.4 


108.8 


0.2 


11.6 


129.1 


111.0 


0.1 


13.5 


126.1 


113.3 


0.0 


Note. 


— The 


circular 


velocities 



are derived under an assumption of a 
variable stellar mass-to-light ratio, as 
described in the text. The CO rota- 
tion model provides a reasonably good 
(though fit by eye) match to the ob- 
served position-velocity diagram. In 
both cases the gas surface density is 
taken to be zero interior to 1.5". 

ferred circular velocity curve. However, at 6.5" (700 pc) 
its rotation velocity is only 70% of the circular velocity 
curve, which means that the enclosed mass is a factor of 
two smaller than that suggested by the circular velocity. 
The discrepancy is even worse at smaller radii, where the 
circular velocity is nearly a factor of three higher than the 
rotation velocities in the molecular gas and the implied 
mass difference is a factor of 10. 

At present the source of the discrepancy between the 
observed CO velocities and the inferred circular velocity 
curve is not understood. The unusually strong peak in 
the circular velocity curve is driven by the very bright 
blue nucleus in optical images and by the high cen- 
tral stellar velocity dispersion. Therefore, one possibil- 
ity is that the measured stellar mean velocity and/or 
velocity dispersion in the center of NGC 3032 are sig- 
nificantly overestimated. However, the velocity dis- 
persion measurements i n the central regions made by 
iMcDermid et al.l (|2006a[ ) with the OASIS instrument are 
consistent with the SAURON data. Another striking 
possibility is that the molecular gas in NGC 3032 may 
not be rotating at the circular velocity. We argued above 
that the gas had to have been acquired from some ex- 
ternal source, and perhaps it is still undergoing strong 
dynamical evolution as it settles into equilibrium. Yet if 
its velocities are < 70% of the circular velocity its orbits 
would be highly noncircular, which seems inconsistent 
with the regular, symmetric appearance of the CO ve- 
locity field (Figure [1]) . Neither of these explanations is 
completely satisfactory. 

In short, the CO velocities in NGC 4459 and NGC 4526 
are in excellent agreement with the circular velocities in- 
ferred for those two galaxies. These matches indicate 
that the stellar and CO kinematics are well understood. 
In NGC 3032 and NGC 4150 at radii > 10", the observed 
CO velocities are nearly equal to the mean stellar veloc- 
ities and 10% to 20% smaller than the inferred circular 
velocities. However, the CO velocities of NGC 3032 are 
conspicuously low when compared to its inferred circu- 
lar velocity within 6" of the nucleus, and it would be 
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valuable to carry out this kind of a comparison in other 
galaxies with counterrotating molecular gas in order to 
ascertain whether NGC 3032 is unusual. 

6. MOLECULAR AND IONIZED GAS VELOCITIES 

In normal spiral galaxies, the bulk of the ionized 
gas emission traces HII regions recently formed out of 
the molecular gas and sharing the disky kinematics of 
the molecular gas. The effect is that in a major axis 
position-velocity diagram the ionized gas tends to trace 
the ridgeline of the molecular and atomic gas or to be 
displaced towards the high velocity envelope (due pri- 
marily to resolution eff e cts). Exa mples can be seen i n 
iKregel fc van der Kruitj (|2004D and lYoung et all (|2006f ). 
Thus, comparisons between the molecular and ionized 
kinematics can help to elucidate the role of star forma- 
tion activity (or the lack thereof) in producing the ion- 
ized gas. 

Excluding the case of NGC 4150, which is discussed 
in S 34.3[ there appears to be a systematic difference be- 
tween ionized gas kinematics and molecular kinematics 
of these lenticulars. The rotation speeds of the ionized 
gas are ~ 20% smaller than those of the CO. There is also 
a slight tendency for the [O III] velocities to be smaller 
than H/3 velocities; this latter feature is most obvious at 
radii > 5" in NGC 4526 (Figure HOI) ■ Lower velocities 
in the ionized gas than in the molecular gas is the oppo- 
site situation to what is usually observed in the disks of 
spirals and is also in the opposite sense to the effects of 
beam smearing on the CO velocities. These data suggest 
that there may be a kinematic component of the ionized 
gas which is not dynamically cold and not related to star 
formation. They also imply that if the ionized gas in 
early-type galaxies is to be used for a measurement of the 
galaxy's circular velocity, an asymmetric drift correction 
maybe significant (e.g. Crctton, Rix, & dc Zccuw 2000). 
By way of a caveat, though, we note that these lentic- 
ulars are dusty and NGC 4526 in particular has a high 
inclination (79°), so the observed CO and ionized gas 
velocities may be biased by projection and optical depth 
effects. The details of the relationships between cold and 
warm ionized phases of the interstellar medium arc not 
yet clear, and they deserve closer scrutiny for the insights 
they may give into the evolution of the ISM in early-type 
galaxies. 

7. SUMMARY 

We present resolved images of the CO emission in 
the four lenticular galaxies NGC 3032, NGC 4150, 
NGC 4459, and NGC 4526. These are some of the 
most CO-rich galaxies in the SAURON survey of early 
type galaxies, so they are prime targets for investigations 
which use cold gas to trace the interaction/merger his- 
tory of early type galaxies and also to document morpho- 
logical change through star formation and disk growth. 
Their inferred H2 masses are in the range 5 x 10^ to 5 x 10^ 
M0. The molecular gas is located in kpc-scale disks (in 
excellent agreement with the distribution of dust visible 
in broadband optical images), the smallest being NGC 
4150 with a radius of 500 pc and the largest being NGC 
3032 with a radius of 1.5 kpc. Average molecular surface 
densities (including helium) are 100 to 200 M0 pc^^. 

In three of the four galaxies (NGC 3032, NGC 4459, 
and NGC 4526) the molecular gas is distributed in disks 



which show regular rotation and little sign of recent dis- 
turbance. The kinematic major axes are well aligned 
with the host galaxies' photometric and (stellar) kine- 
matic major axes, suggesting that the gas has settled into 
the equatorial plane of nearly axisymmetric potentials. 
The velocity field of NGC 4150 shows a kinematic major 
axis in rough agreement with the galaxy's optical major 
axis, but better spatial resolution will be necessary in 
order to assess the regularity of the CO kinematics. Fur- 
thermore, in NGC 3032 the molecular gas's kinematic 
position angle is 180° offset from the stellar kinematic 
position angle; this dramatic counterrotation indicates 
that the molecular gas was acquired from an external 
source or perhaps is leftover from a major merger, but it 
cannot have been produced through internal stellar mass 
loss. The sense of the CO rotation in NGC 3032 is con- 
sistent with that of the young kinematically decoupled 
core, however, which suggests that this is an example of 
a stellar substructure forming through dissipational pro- 
cesses. 

In two cases (NGC 4459 and NGC 4526) the CO 
kinematics provide powerful, independent confirma- 
tion of the dyn a mical mass-to-light ratios inferred by 
ICappellari et al.l ()2006[ ). The mass-to-light ratios were 
derived from the SAURON stellar kinematic data via 
two-integral Jeans and three-integral Schwarzschild dy- 
namical models and were used to infer circular velocity 
curves. Simple tilted ring models are presented in which 
the molecular gas rotates at the circular velocity (or, in 
NGC 4526, 95% of the circular velocity, a decrease which 
is roughly consistent with the uncertainty in the mass- 
to-light ratio). Comparison of the model and observed 
major axis position-velocity diagrams indicates that the 
circular velocity inferred from stellar kinematics alone is 
consistent with the behavior of the molecular gas. The 
agreement is significant because it is rare to be able to 
make this kind of an independent check of the stellar 
dynamical analysis. 

Several puzzling results of this work deserve further 
study; these are disagreements between the gaseous and 
the stellar kinematics at sub-kpc scales. For example, 
in NGC 4150 the CO velocities appear to be inconsistent 
with those of the young (counterrotating) stellar core. In 
NGC 3032 the observed CO velocities are significantly 
lower than the circular velocity curve (amounting to an 
implied factor of two less mass interior to 700 pc and 
even more at smaller radii). At present it is not known 
whether the stellar kinematic data give an overestimated 
circular velocity for this galaxy or whether the recently 
accreted molecular gas is not following circular orbits. 
Higher resolution CO observations will also be neces- 
sary to check whether NGC 4150 shows an inconsistency 
between CO velocities and circular velocities. Finally, 
the relationships between molecular gas and ionized gas 
should contribute valuable insights into star formation 
and the evolution of the ISM in these lenticulars. 
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